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The  coherence  of  an  optical  beam  having  multiple  degrees  of  freedom  (DoFs)  is  described  by  a 
coherency  matrix  G  spanning  these  DoFs.  This  optical  coherency  matrix  has  not  been  measured 
in  its  entirety  to  date— even  in  the  simplest  case  of  two  binary  DoFs  where  G  is  a  4  x  4  matrix. 

We  establish  a  methodical  yet  versatile  approach— optical  coherency  matrix  tomography— for 
reconstructing  G  that  exploits  the  analogy  between  this  problem  in  classical  optics  and  that  of 
tomographically  reconstructing  the  density  matrix  associated  with  multipartite  quantum  states  in 
quantum  information  science.  Here  G  is  reconstructed  from  a  minimal  set  of  linearly  independent 
measurements,  each  a  cascade  of  projective  measurements  for  each  DoF.  We  report  the  first 
experimental  measurements  of  the  4X4  coherency  matrix  G  associated  with  an  electromagnetic 
beam  in  which  polarization  and  a  spatial  DoF  are  relevant,  ranging  from  the  traditional  two-point 
Young's  double  slit  to  spatial  parity  and  orbital  angular  momentum  modes. 

The  statistical  fluctuations  of  light  in  space  and  time  may  be  characterized  by  a  hierarchy  of  correlation 
functions  for  electromagnetic  field  components1,2.  These  functions,  not  the  optical  fields  themselves,  pro¬ 
vide  a  description  of  light  in  terms  of  observable  quantities3.  The  theory  of  optical  coherence  investigates 
the  properties  of  these  correlation  functions  pertaining  to  the  temporal,  spatial,  and  polarization  degrees 
of  freedom  (DoFs).  When  these  DoFs  are  uncoupled  (or  uncorrelated),  simple  measures  of  coherence  for 
each  DoF  suffice,  such  as  coherence  time,  coherence  area,  and  degree  of  polarization4.  However,  when 
the  DoFs  are  coupled ,  such  measures  lose  their  utility  and  more  sophisticated  approaches  are  required, 
such  as  the  mutual  coherence  function5,  the  beam  coherence-polarization  matrix6-8,  or  the  4x4  field 
correlation  matrix  for  a  pair  of  points  in  an  electromagnetic  field9-11. 

While  the  importance  of  coupling  between  DoFs  was  recognized  decades  ago,  as  in  Mandels  sem¬ 
inal  work  on  optical  cross-spectral  purity  (the  absence  of  spatial-spectral  coupling)12,  recent  advances 
have  led  to  a  host  of  scenarios  wherein  such  coupling  is  critical.  For  example,  vector  beams  correlate 
polarization  with  spatial  position13,  scattering  from  complex  photonic  structures  and  devices  may  cou¬ 
ple  the  relevant  field  DoFs14,15,  and  reliance  on  multimode  optical  fibers  for  spatial  multiplexing  is 
reviving  interest  in  joint  polarization-spatial-mode  characterization16.  In  exploring  these  settings,  it  has 
recently  proven  fruitful  to  adopt  the  Hilbert- space  formulation  used  in  quantum  mechanics  to  the  needs 
of  classical  coherence  theory10,11 — an  approach  that  has  early  prescient  antecedents17,18.  In  the  context 
of  coupling  between  multiple  DoFs,  such  a  treatment  necessitates  introducing  the  notion  of  classical 
entanglement’10,19-25.  In  quantum  mechanics,  states  associated  with  bipartite  systems  that  do  not  sep¬ 
arate  into  products  of  states  belonging  to  the  Hilbert  space  of  each  particle  are  said  to  be  entangled26. 
As  a  consequence  of  the  mathematical  similarity  between  the  Hilbert  spaces  of  multi-partite  quantum 
states  and  multi-DoF  classical  optical  fields,  a  corresponding  concept  of  classical  entanglement  indi¬ 
cates  the  non-separability  of  the  beam  into  uncoupled  DoFs.  After  the  initial  suggestion  by  Spreeuw19, 
a  substantial  body  of  work  has  accumulated  in  the  past  five  years  in  which  classical  entanglement  is 
exploited  in  solving  long-standing  problems  in  polarization  optics27-29,  delineating  the  contributions  of 
non-separability  and  intrinsic  randomness  to  the  coherence  of  an  optical  beam10,30,  introducing  new 
metrology  schemes31,  and  implementing  classical  analogs  of  quantum  information  processing  protocols, 
such  as  teleportation32,33,  and  super-dense  coding,  etc. 
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A  fundamental  capability  that  has  remained  elusive  in  classical  optics  is  the  complete  identification 
of  the  coherence  function  for  a  beam  with  coupled  DoFs.  In  quantum  mechanics,  the  task  of  measuring 
all  the  elements  of  a  density  matrix  is  known  as  quantum  state  tomography’34,35.  The  corresponding 
procedure  for  multi-DoF  beams  in  classical  optics  has  been  studied  theoretically11,  but  has  not  been 
demonstrated  experimentally  heretofore.  Even  in  the  simplest  case  of  two  binary  DoFs6  (e.g.,  polariza¬ 
tion,  a  bimodal  waveguide36,37,  two  coupled  single-mode  waveguides38,39,  spatial-parity  modes40-44,  etc.), 
the  associated  4x4  coherency  matrix  G,  which  is  a  complete  representation  of  second-order  coher¬ 
ence10,11,  has  not  been  measured  in  its  entirety  to  date. 


Results 


In  this  Article,  we  present  a  methodical  approach — optical  coherency  matrix  tomography 
(OCmT) — for  measuring  the  complex  elements  of  4  x  4  coherency  matrices  G  by  appropriating  the 
quantum-state-tomography  strategy.  To  demonstrate  the  universality  of  our  approach,  we  implement 
it  with  coherent  and  partially  coherent  fields  having  coupled  or  uncoupled  DoFs  in  three  distinct  set¬ 
tings  involving  pairs  of  points9-11,  spatial-parity  modes40-44,  and  orbital  angular  momentum  (OAM) 
modes45 — each  together  with  polarization.  We  identify  the  minimal  set  of  linearly  independent,  joint 
spatial-polarization  projective  measurements  that  enable  a  unique  reconstruction  of  G.  Since  G  is  a  com¬ 
plete  representation  of  the  field,  its  reconstruction  obviates  the  need  to  measure  directly  any  coherence 
descriptors  (all  of  which  are  scalar  functions  of  the  complex  elements  of  G)  and,  moreover,  allows  for 
unambiguous  identification  of  classical  entanglement. 

The  coherence  of  an  optical  beam  having  a  single  binary  DoF  is  represented  by  a  2  x  2  Hermitian 


coherency  matrix  G  = 
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to  one  level  of  the  DoF.  For  example,  polarization  is  represented  by  the  coherency  matrix4 
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where  Gy,  =  (?) E  y(j)  Ej  (?)  is  the  horizontal  (H)  or  vertical  (V)  field  component  at  a  point  ?,  {Sf} 

are  the  Stokes  parameters,  and  {d/}  are  the  Pauli  matrices 


Polarization  coherence  is  quantified  by  the  degree  of  polarization  D  —  +  (S|)2  +  (S3P)2 

— with  Dp  =  1  and  Dp  =  0  corresponding  to  purely  polarized  and  completely  unpolarized  light,  respec¬ 
tively.  The  Stokes  parameters  are  evaluated  via  four  projective  polarization  measurements  [Fig.  1(a)]:  If, 
If,  and  If  correspond  to  the  H,  diagonal  (D),  and  right-hand-circular  (R)  polarization  components, 
respectively,  in  addition  to  the  total  power  If  ref.  46;  in  which  case  Sf  =  2 If  —  If.  The  same  formalism 
may  be  applied  to  other  binary  DoFs  [Fig.  1(b)]:  (i)  the  scalar  field  at  two  points  ?a  and  ?b,  £a  and  Eb; 
(ii)  the  spatial-parity  even  e’  and  odd  o’  modes  of  a  scalar  field  Ee  and  E0  refs.  40-44;  or  (iii)  a  pair  of 
OAM  modes,  e.g.,  E0  and  El  corresponding  to  OAM  £  =  0  and  1,  respectively45,47. 

When  two  binary  DoFs  of  the  field  are  relevant,  e.g.,  the  first  is  polarization  p’  and  the  second  is  a 
spatial  s’  DoF  with  modes  identified  as  a’  and  £b’  [Fig.  1(c)]— the  corresponding  coherency  matrix  G  is 
now  4x4  refs.  10,  11, 


G  = 


GHa,Ha 

GHa,Hb 

GHa,Va 

GHa,Vb 

GHb,Ha 

GHb,Hb 

GHb,Va 

GHb,vb 

1  A  c  ~p 

Gva,Ha 

Gya,Hb 

Gya,Va 

Gya,Vb 

=  7  E  Sl,mtf 

4  /,m= 0 

®  Pm, 

Gyb,Ha 

Gyb,Hb 

Gyb,Va 

Gyb,Vby 

' )’  Gjk,j'k' 

=  Gfk'. 

is 

a  field  component  with  j 

.  C  - 

=  H,  V 

(3) 


the  two -DoF  Stokes  parameters,  and  {dj3}  and  {d^}  are  the  Pauli  matrices  on  the  polarization-  and 
spatial-DoF  Hilbert  subspaces,  respectively11.  In  determining  the  coherence  descriptors  of  each  DoF 
independently  of  the  other,  one  first  traces  over  the  other  DoF  to  obtain  a  reduced  coherency  matrix10. 
The  reduced  polarization  coherency  matrix  Gp,  obtained  by  tracing  over  the  spatial  DoF  in  G,  is  given 
by 


GHa,Ha  +  GHo?Ho  GHaVa  +  GHb,Vb 
Gya,Ha  +  GVb,Hb  GVaVa  +  G  vbjVb 


(4) 
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Figure  1.  Measurement  scheme  for  optical  coherency  matrix  tomography  (OCmT).  (a)  Four 
projections  If,l  =  0,  ...,  3  to  obtain  the  Stokes  parameters  {St}  and  reconstruct  Gp  for  the  polarization 
DoF.  (b)  Similarly,  four  projections  I^,m  =  0,  ...,  3  to  obtain  the  Stokes  parameters  {Sm}  and  reconstruct 
Gs  for  a  binary  spatial  DoF.  (c)  OCmT  enables  the  reconstruction  of  G  for  the  two  binary  DoFs  in 
(a,b)  considered  simultaneously  via  16  joint  polarization- spatial  measurements,  each  of  which  consists  of  a 
cascade  of  two  projections— one  from  (a)  and  the  other  from  (b). 


while  the  reduced  coherency  matrix  for  the  spatial  DoF  Gs,  obtained  after  tracing  over  polarization  in  G 
(see  ref.  11),  is  given  by 


^Ha,Ha  +  ^Va,Va  ^Ha,Hb  +  ^Va,Vb 
^Hb,Ha  +  ^Vb,Va  ^Hb,Hb  +  ^Vb,Vb^ 


(5) 


The  elements  of  the  reduced  coherency  matrices  are  measured  by  a  system  sensitive  to  one  DoF,  but 
not  to  the  other.  When  the  two  DoFs  are  uncoupled,  then  G  =  Gp  0  Gs,  otherwise  the  elements  of  Gp 
and  Gs  lack  information  about  the  correlations  between  the  two  DoFs  that  is  contained  in  G.  Such  cor¬ 
relations  are  only  measurable  by  a  system  that  is  sensitive  to  both  DoFs  via  joint  polarization-spatial 
measurements. 

We  pose  the  following  question:  what  are  the  necessary  and  sufficient  measurements  to  reconstruct 
an  arbitrary  G  for  two  binary  DoFs?  This  question  was  solved  by  Wootters48  in  the  context  of  recon¬ 
structing  the  density  matrix  p  for  a  bipartite  quantum  system.  He  showed  that  the  measurements  carried 
out  on  each  subsystem  to  reconstruct  its  reduced  density  matrix  are  sufficient  to  reconstruct  p  when 
carried  out  jointly — a  methodology  known  as  quantum  state  tomography34,35.  In  our  context  of  a  classical 
optical  beam  having  two  binary  DoFs,  the  analogy  with  the  quantum  setting  allows  us  to  exploit  the  same 
strategy.  Regardless  of  the  specific  form  of  G,  the  necessary  measurements  to  carry  out  OCmT  and 
reconstruct  G  [Fig.  1(c)]  are  those  used  to  reconstruct  the  reduced  coherency  matrices  [Fig.  l(a,b)] 
carried  out  in  cascades  of  pairs  of  projections — one  for  polarization  and  the  other  for  the  spatial  DoF. 
Each  measurement  yields  a  real  number  IXm  (projection  l  for  polarization  and  m  for  the  spatial  DoF) 
corresponding  to  the  projection  of  a  tomographic  slice  through  G.  The  16  combinations  of 
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Figure  2.  Polarization  considered  with  spatial  location  in  a  Young’s  double-slit  configuration. 

(a)  Experimental  setups  for  preparing  polarization  and  position-coupled  beams  Gx  through  G6;  D:  diffuser 
that  randomizes  the  beam  spatially;  PS:  polarization  scrambler  that  randomizes  the  beam  polarization; 

HWP:  half-wave  plate,  (b)  Experimental  setup  delineating  the  stages  of  polarization  analysis  followed  by 
spatial  analysis.  HWP:  half-wave  plate;  QWP:  quarter-wave  plate;  PBS:  polarizing  beam  splitter;  BC:  50:50 
beam  coupler;  PD:  phase  delay  element  that  introduces  a  phase  shift  ip.  (c)  Illustration  depicting  the 
equivalence  between  the  phase  shift  ip  introduced  by  the  phase  element  PD,  and  the  lateral  displacement  x 
on  a  screen  upon  which  the  far-held  intensity  pattern  is  projected,  (d)  Spatial  profile  measurements  obtained 
by  a  CCD  camera  illustrating  the  spatial  projective  measurements  for  the  H  polarization  projection;  similarly 
for  the  V,  D,  and  R  projections. 


polarization-spatial  measurements  are  inverted  to  obtain  the  two-DoF  Stokes  parameters, 
Sim  =  4 Ilm  —  2 Il0  —  2 10m  +  700,  and  hence  reconstruct  G  (see  ref.  11  for  details). 

We  have  performed  a  series  of  experiments  implementing  the  OCmT  scheme  described  above  using 
quasi-monochromatic  beams  having  two  binary  DoFs:  polarization  and  a  spatial  DoF.  We  have  measured 
the  4x4  coherency  matrix  G  for  six  different  beams  corresponding  to  distinct  states  of  light  having  the 
following  properties: 

Gp  the  polarization  and  spatial  DoFs  are  separable  and  both  are  coherent. 

G2:  the  polarization  DoF  is  coherent  while  the  spatial  DoF  lacks  coherence. 

G3:  both  the  polarization  and  spatial  DoFs  lack  coherence. 

G4:  the  polarization  and  spatial  DoFs  are  classically  entangled. 

G5:  the  polarization  and  spatial  DoFs  are  classically  correlated. 

G6:  this  beam  is  a  mixture  of  the  separable-coherent  beam  Gx  and  the  classically  entangled  beam  G4. 

We  use  the  sequence  of  polarization  projections  described  earlier  and  present  below  the  spatial  pro¬ 
jections  following  the  H  projection  (similar  spatial  projections  are  carried  out  following  the  V,  D  and  R 
polarization  projections). 

Polarization  with  Spatial  Position.  The  first  realization  of  the  spatial  DoF  is  the  traditional  two 
points,  as  in  the  Youngs  double  slit  experiment.  The  polarization  and  position-coupled  beam  is  prepared 
in  one  of  six  states  Gx  through  G6;  Fig.  2(a)  (see  Supplementary  information  for  details).  OCmT  for  a 
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Figure  3.  Polarization  considered  with  spatial  parity  modes,  (a)  Experimental  setups  for  preparing 
polarization  and  spatial-parity-coupled  beams  Gx  through  G6;  HWP:  half-wave  plate;  SLM:  spatial  light 
modulator;  PS:  polarization  scrambler  that  randomizes  the  beam  polarization,  (b)  Experimental  setup 
delineating  the  stages  of  polarization  analysis  followed  by  spatial-parity  analysis.  HWP:  half-wave  plate; 
QWP:  quarter-wave  plate;  PBS:  polarizing  beam  splitter;  SLM:  spatial  light  modulator;  BS:  beam  splitter;  SF: 
spatial  flipper  that  flips  the  sign  of  the  odd  ‘o’  spatial-parity  mode,  leaving  the  even  e’  mode  intact. 

(c)  Corresponding  spatial  profile  measurements  obtained  by  a  CCD  camera  illustrating  the  spatial  projective 
measurements  for  the  H  polarization  projection;  similarly  for  the  V,  D,  and  R  projections. 


such  a  beam  comprises  of  the  polarization  analysis  followed  by  the  spatial  analysis;  Fig.  2(b).  The  spatial 
analysis  may  alternatively  be  carried  out  by  extracting  specific  intensity  points  from  the  far-held  intensity 
patterns  for  only  two  values  of  displacement  x  on  a  screen  or  an  array  of  detectors;  Fig.  2(c).  The  four 
spatial  projections  are  obtained  by  measuring  the  following:  (1)  the  total  power  from  both  points  ho  — 
at  x  =  0;  (2)  the  power  from  point  a ’  In  =  JHa  at  x=  0;  (3)  the  power  in  the  far-held  interference  pattern 
I \2  —  f H(a+b)  at  *=  0;  and  (4)  the  power  at  the  value  of  x  corresponding  to  a  -  phase  shift  I13  =  fH(a+/b); 
see  Fig.  2(d).  It  is  important  to  note  that  the  visibility  of  fringes  is  not  the  parameter  sought  here  to 
characterize  the  spatial  coherence  at  ‘a  and  £b’;  instead  the  four  points  identihed  in  Fig.  2(d),  together 
with  the  set  of  points  obtained  for  the  V,  D,  and  R  polarization  projections,  reveal  the  complete  picture 
even  when  polarization  and  the  spatial  DoFs  are  classically  entangled. 

Polarization  with  Spatial  Parity.  The  second  spatial-DoF  realization  makes  use  of  one-dimensional 
even  e’  and  odd  o’  spatial-parity  modes  with  respect  to  x=  0.  The  polarization  and  spatial  parity- coupled 
beam  is  prepared  in  one  of  six  states  Gx  through  G6;  Fig.  3(a)  (see  Supplementary  information  for  details). 
OCmT  for  a  such  a  beam  comprises  of  the  polarization  analysis  followed  by  the  spatial-parity  analysis; 
Fig.  3(b).  The  four  spatial  projections  are  obtained  by  measuring  the  power  (integrated  over  the  shaded 
areas  in  Fig.  3(c))  in  the  following  settings:  (1)  the  total  power  ho  —  Ih  of  the  beam;  (2)  the  power  of  the 
even  component  hi  —  obtained  from  a  modified  Mach-Zehnder  interferometer  that  separates  the 
beam  into  the  different  spatial-parity  components41;  (3)  the  power  Iu  =  JH(e+0)  after  blocking  the 
half-plane  x<  0,  corresponding  to  a  projection  onto  the  e  +  o  component;  and  (4)  a  projection  onto  the 
e  +  io  component  obtained  from  the  power  J13  =  IH(e+fo)  even  component  measured  after  first 
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Figure  4.  Polarization  considered  with  OAM  modes,  (a)  Experimental  setups  for  preparing  polarization 
and  OAM-coupled  beams  Gx  through  G6;  HWP:  half-wave  plate;  SLM:  spatial  light  modulator;  PS: 
polarization  scrambler  that  randomizes  the  beam  polarization,  (b)  Experimental  setup  delineating  the  stages 
of  polarization  analysis  followed  by  OAM-mode  analysis.  HWP:  half-wave  plate;  QWP:  quarter-wave  plate; 
PBS:  polarizing  beam  splitter;  SLM:  spatial  light  modulator;  FC:  fiber  coupler;  SMF:  single-mode  fiber, 
(c)  Corresponding  spatial  profile  measurements  obtained  by  a  power-meter  illustrating  the  spatial  projective 
measurements  for  the  H  polarization  projection;  similarly  for  the  V,  D,  and  R  projections.  Spatial 
measurements  are  obtained  by  dislocating  the  phase  singularity  (e2<^,  0  <  0  <  27t)  relative  to  the  beam  Ax 
along  x  and  Ay  along  y.  IH  is  obtained  by  adding  the  intensities  Ie  at  Axmax  a  (a  is  the  beam  width)  and 
I£  at  Ax=  0,  I£  +i  is  obtained  at  Axmid,  and  I£  +ii  is  obtained  at  Aymid  (Axmid  and  Aymid  are  calibrated 
using  a  Gaussian  beam). 


introducing  a  phase-step  -  between  the  two  plane  halves  x  <  0  and  x  >  0  implemented  by  a  spatial  light 
modulator  (SLM);  see  Fig.  3(c).  This  phase  modulation  was  shown  in  ref.  40-44  to  produce  a  rotation 
on  a  major  circle  on  a  Poincare  sphere  having  the  e  and  o  modes  as  antipodes. 

Polarization  with  OAM.  The  third  realization  exploits  two  low-order  OAM  modes  £  =  0  (£0)  and 
£  =  1  (£l).  The  polarization  and  OAM-coupled  beam  is  prepared  in  one  of  six  states  Gl  through  G6; 
Fig.  4(a)  (see  Supplementary  information  for  details).  OCmT  for  a  such  a  beam  comprises  of  the  polar¬ 
ization  analysis  followed  by  the  OAM  mode  analysis;  Fig.  4(b). The  four  spatial  projections  are  obtained 
by  measuring  the  following:  (1)  the  total  power  (2)  the  power  of  the  £0  component  710  =  IHi 

obtained  using  a  spatial  filter  (a  lens  focusing  to  a  single-mode  fiber);  (3)  the  power  of  the  £0  +  £x  com¬ 
ponent  using  a  phase  vortex  with  the  dislocation  displaced  laterally  with  respect  to  the  beam  imple¬ 
mented  using  an  SLM  Il2  =  IH(T0+q)  re^-  49;  and  (4)  the  power  of  the  £0  +  i£x  component  using  the  same 
phase  vortex  but  with  the  dislocation  displaced  vertically  J13  =  IH(4+/£1);  see  Fig-  4(c). 

Measurements.  We  have  measured  the  complex  elements  of  G  for  six  different  classes  of  beams 
comprising  those  with  separable  DoFs  (both  coherent,  both  incoherent,  or  in  a  hybrid  coherent/ incoher¬ 
ent  configuration),  non-separable  DoFs  (classically  entangled  or  classically  correlated),  and  mixtures  of 


SCIENTIFIC  REPORTS  |  5:15333  |  DOI:  io.i038/srepi5333 


6 


www.nature.com/scientificreports/ 


Figure  5.  Measurements  of  G  for  beams  having  separable  DoFs,  (a)  Gp  both  DoFs  are  coherent;  (b)  G2: 
polarization  is  coherent  but  the  spatial  DoF  is  incoherent;  and  (c)  G3:  both  DoFs  are  incoherent.  In  each 
case,  the  reduced  coherency  matrices  Gp  and  Gs  are  also  depicted.  The  imaginary  components  are  negligible, 
and  are  not  shown. 


beams  from  the  two  separable  and  non-separable  classes  (see  Supplementary  information  for  the  com¬ 
plete  results).  In  each  experiment,  the  prepared  beam  passes  first  through  polarization  then  spatial-DoF 
analysis  stages  (the  order  may  be  reversed  without  changing  the  outcome).  In  each  of  these  realizations, 
permutations  of  the  four  polarization  projection  settings  combined  with  the  four  spatial  projection  set¬ 
tings  yield  16  measurements  for  OCmT,  which  are  used  to  reconstruct  G.  We  make  use  of  a 
maximal-likelihood  algorithm  that  exploits  the  constraints  set  by  the  trace,  hermiticity,  and 
semi-positive-definiteness  of  G50.  We  portray  the  real  and  imaginary  components  of  G  using  the  standard 
visualization  from  quantum  state  tomography.  In  each  plot  we  provide  the  coherence  descriptor  for  the 
polarization  Dp  and  spatial  DoF  Ds  obtained  from  their  reduced  coherency  matrices,  in  addition  to  the 
linear  entropy  DG  =  ^  |Tr(G2)  —  ^|,  which  serves  as  a  measure  of  the  overall  beam  coherence,  where 
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Figure  6.  Measurements  of  G  for  beams  having  non-separable  DoFs,  (a)  G4:  classically  entangled  beam; 
and  (b)  G5:  classically  correlated  beam.  The  imaginary  components  in  both  cases  are  negligible,  and  are  not 
shown,  (c)  Measurement  of  the  real  and  imaginary  parts  of  G6:  mixture  of  beams  Gx  and  G4.  In  each  case, 
the  reduced  coherency  matrices  Gp  and  Gs  are  also  depicted. 


Dg  ranges  from  0  (complete  incoherence  in  all  DoFs)  to  1  (a  coherent  beam  with  no  statistical  fluctua¬ 


tions)51.  Finally,  we  provide  the  fidelity  F  = 


TrjVVrGVr  | 


as  a  measure  of  the  robustness  of  the 


reconstruction  process  via  OCmT52,  where  T  is  the  theoretical  matrix  and  G  is  the  measured  matrix. 


Beams  with  separable  DoFs.  We  present  in  Fig.  5  three  examples  of  beams  having  separable  DoFs, 
G  =  Gp  0  Gs.  First,  both  polarization  and  spatial  DoFs  are  coherent,  Gx  =  (o  o)  ®  i(i  ;)[*§•  *a)]- 
Second,  a  hybrid  beam  in  which  polarization  is  pure  (along  D)  but  the  beam  is  spatially  incoherent, 

°2  =  Ki  l)  0  Ho  ° )  [pis- Ihird> a  comPletely incoherent  beamG3  =  i)  ®  i(o  i) 

[Fig.  5(c)].  In  all  three  cases,  the  separability  of  the  two  DoFs  is  readily  detected  by  visual  inspection  of  G 
and  confirmed  by  taking  the  direct  product  of  the  reduced  coherency  matrices. 

Beams  with  non-separable  DoFs.  We  next  present  two  fundamentally  distinct  classes  of  beams  with 
non-separable  G  in  Fig.  6(a,b).  First,  OCmT  of  a  classically  entangled  beam  G4  is  shown  in  Fig.  6(a), 
wherein  the  beam  is  fully  coherent ,  and  yet  the  measures  extracted  from  reduced  coherency  matrices 
indicate  complete  incoherence.  In  such  a  beam,  the  polarization  and  spatial  modes  occur  in  pairs— e.g., 
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H  with  ‘a  and  Y  with  £b’  (but  never  H  with  £b’  or  Y  with  £a).  In  the  traditional  view  of  the  double  slit 
experiment,  such  coupling  will  produce  no  interference  fringes,  and  the  lack  of  visibility  may  be  inter¬ 
preted  as  the  absence  of  spatial  coherence,  despite  the  beam  being  perfectly  spatially  coherent.  This 
coupling  between  the  DoFs  is  in  fact  encoded  in  the  non-zero  off-diagonal  elements  of  G  revealed  once 
it  is  reconstructed  through  OCmT,  but  cannot  be  obtained  from  Gp  or  Gs.  Second,  a  classically  correlated 
beam  G5  is  shown  in  Fig.  6(b)  in  which  the  same  coupling  between  polarization  and  spatial  modes  occurs 
as  in  the  previous  example,  except  the  different  combinations  are  incoherently  mixed  and  not  linearly 
superposed.  The  partial  global  coherence — despite  the  complete  lack  of  coherence  for  each  DoF — is  clear 
from  the  fact  that  not  all  the  diagonal  elements  of  G5  are  equal  as  is  the  case  in  G3. 

Mixture  of  beams  with  separable  and  non-separable  DoFs.  Finally,  in  Fig.  6(c)  we  depict  G6  correspond¬ 
ing  to  a  beam  formed  by  statistically  mixing  the  separable-coherent  beam  Gx  and  the  classically  entan¬ 
gled  beam  G4.  The  measurement  of  G6  indicates  that  part  of  the  apparent  incoherence  in  this  beam  stems 
from  the  intrinsic  randomness  in  the  individual  DoFs,  and  part  of  it  from  the  correlation,  or  classical 
entanglement,  between  the  two  DoFs. 

Discussion 

The  reconstruction  of  G  allows  for  the  unambiguous  and  complete  mathematical  expression  of  fields 
that  are  coherent,  partially  coherent,  or  incoherent,  in  either,  or  both,  DoFs  of  an  optical  beam  with  two 
binary  DoFs.  The  usefulness  of  this  technique  becomes  specially  apparent  in  cases  where  the  DoFs  are 
coupled  or  non-separable,  and  the  traditional  scalar  measures  of  coherence  provide  a  conflicting  and 
fallacious  account  of  beam  coherence.  The  apparent  absence  of  coherence  in  any  DoF  may  be  the  result 
of  intrinsic  randomness  due  to  statistical  fluctuations,  or  due  to  the  coupling  or  non-separability  with 
another  DoF.  In  the  latter  case,  the  measurement  of  G  also  provides  the  way  for  implementing  unitary 
transformations  required  to  undo  such  coupling,  and  restore  coherence  in  the  DoFs.  The  application  of 
our  work  can  be  easily  seen  in  the  myriad  applications  of  coherence  under  conditions  of  coupled  DoFs, 
particularly  those  involving  localized  vector  beams,  sub -diffraction  imaging,  nanophotonics,  and  prop¬ 
agation  through  disordered  media.  Measurement  of  G,  before  and  after  transmission  though  a  system 
that  couples  various  DoFs,  will  help  determine  the  characteristics  of  the  system.  This  technique  may 
hence  find  important  applications  in  crystallography,  atmospheric  optics,  and  systems  involving  photonic 
crystals  or  anisotropic  scatters,  etc. 

In  summary,  we  have  experimentally  demonstrated  for  the  first  time  a  methodical,  yet  versatile, 
approach  to  reconstructing  the  4x4  coherency  matrix  G  of  an  optical  beam  having  two  binary  DoFs, 
which  we  call  optical  coherency  matrix  tomography.  We  have  explored  three  different  physical  realizations 
in  which  we  combine  polarization  with  spatial  position,  spatial  parity,  or  orbital  angular  momentum 
modes.  By  exploiting  the  mathematical  similarity  with  quantum  state  tomography  of  two  photon  states, 
we  determine  the  minimal  set  of  measurements  required  to  reconstruct  G.  Although  we  have  conducted 
the  experiments  for  a  beam  with  two  binary  DoFs,  this  methodology  is  equally  applicable  for  a  higher 
number  of  DoFs  with  m-ary  levels  each. 
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